Cavity QED with magnetically coupled collective spin states 
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We report strong coupling between an ensemble of nitrogen-vacancy center electron spins in 
diamond and a superconducting microwave coplanar waveguide resonator. The characteristic scaling 
of the collective coupling strength with the square root of the number of emitters is observed directly. 
Additionally, we measure hyperfine coupling to 13 C nuclear spins, which is a first step towards a 
nuclear ensemble quantum memory. Using the dispersive shift of the cavity resonance frequency, we 
measure the relaxation time of the NV center at mK temperatures in a non-destructive way. 
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Building a practical quantum information processor 
will profit from a hybridization of different quantum sys- 
tems in order to fulfill the requirements of long stor- 
age times, fast processing speeds, long-distance infor- 
mation transfer and scalability. A variety of hybrid 
systems have been conceived, including cold atoms and 
ions, nanomechanics, and molecules [1—6]. Superconduct- 
ing (SC) qubits [7] are appealing processing devices due 
to the fast processing speed and their customizability. 
Superconducting coplanar waveguide (CPW) resonators 
have been used to reversibly transfer quantum states be- 
tween two qubits and the resonator [8]. Color centers in 
solids on the other hand - most prominently the nega- 
tively charged nitrogen-vacancy defect in diamond (NV) 
- show astonishingly long coherence times even at room 
temperature [9] and provide coherent bridges between 
electron spin resonances (ESR) in the GHz range and 
optical photons suitable for long-distance transfer [10]. 

The coupling between a single spin and the electro- 
magnetic field is typically rather weak. However, when 
writing single excitations into ensembles of N spins, it 
is enhanced by a factor y/N [11]. Provided that excita- 
tions can be coherently written to and retrieved from the 
ensemble by strong coupling to high-Q resonators, quan- 
tum memories can be constructed [12, 13]. Efforts paral- 
lel to this work have succeeded in demonstrating strong 
ensemble-resonator coupling using NVs [14] and other 
impurity spins [15-17]. Schemes such as gradient echoes 
[18] and controlled reversible inhomogeneous broadening 
[19] can be used to store quantum information [20]. 

In this Letter, we employ a strongly coupled system 
consisting of a CPW resonator and an ensemble of NVs to 
directly demonstrate key features of ensemble coupling, 
such as the v^V scaling of the coupling rate. Properties 
of the spin ensemble including the NV spin-lattice relax- 
ation time Ti and ESR zero-field splitting parameters, 
are determined via the cavity resonance shift caused by 
dispersive interactions. We also observe the hyperfine 



interaction between NVs and nearby 13 C nuclear spins. 

The experimental setup is sketched in Fig. 1: A 
4.5 x 2.25 x 0.5 mm 3 high-temperature high-pressure 
diamond is placed on top of a superconducting copla- 
nar waveguide resonator and fixated via a spring-loaded 
metallic pin. Our single crystal diamond sample has two 
polished (100) surfaces and is specified to contain a nitro- 
gen impurity concentration of about 200 ppm. In order 
to form negatively charged NV defect centers, the sample 
was neutron irradiated with a fluence of 5 x 10 17 cm -2 
to create vacancies [21]. Subsequent annealing for four 
hours at 900 °C resulted in an NV concentration of about 
6 ppm which was determined by comparing the fluo- 
rescence intensity of an ensemble to a single NV. Res- 
onator chip and diamond are enclosed in a copper box 
and positioned in the center of two pairs of perpendicular 
Helmholtz coils. This coil configuration allows to apply 
a homogeneous magnetic field B (up to about 25 mT) 
over the extent of the diamond and allows for adjusting 




FIG. 1. (a) (Color online) Schematic of the experimental 
setup. A (001) diamond is positioned in the middle and par- 
allel to the surface of a A/2 CPW resonator. The chip di- 
mensions are 12 x 4 mm, the center conductor is 20 /im wide, 
the gaps are 8.3 /im. A 2-axis Helmholtz coil configuration 
creates a magnetic field in an arbitrary direction ip within the 
(001) plane, tuning the NV centers, (b) Sketch of a diamond 
lattice with the four possible NV center orientations. Due to 
the diamond symmetry, the NV ensembles denoted as I (red) 
and II (green) are tuned equivalently, since these ensembles 
enclose the same angle with the applied magnetic field. 



2 



, (c) ^10-3 




^^h^^^^^^^^^H^bh^^^^H^^^^B 1 1 2 4 £ 

5 10 10 15 20 5 10 15 20 5 10 15 20 5 10 15 20 

B(mT) B(mT) B (mT) B (mT) B (mT) 



FIG. 2. (Color online) Resonator transmission jS^ij 2 as a function of the magnetic field for (a) ip — 45° and (b) tp — 3°. While 
in (a) an avoided crossing of the coupled resonator-spin ensemble system is observed when u) r = ujL , we see in (b) the additional 
splitting of the w 1 ] transition, (c) Comparison of the vacuum Rabi splitting of 2gj /2ir — 18.5 MHz for tp — 45° (B — 8.2 mT) 
(gray line) and the \[2 larger splitting of 2g/2n — 26.3 MHz (black line) for tp — 0° (B = 14.3 mT). Each dataset is fitted by 
a two-peak Lorentzian curve, (d) Transition frequencies uj± and ui± between the ground |0}_r t j7 and the excited states \±}ui 
for different magnetic field angles (tp — 45°, 3°, 0°). The dashed lines indicate the cavity frequency. 



the magnetic field direction in the plane of the resonator 
surface, without changing the resonators frequency [22]. 

To explain our experimental observations, we consider 
the spin Hamiltonian of a single NV center [23] 

H NV = gePeB • S + SDS - g n /3 n B • I + IAS, (1) 

where the first and third term describe the Zeeman inter- 
action of an electronic NV spin (S = 1, g e = 2.0028) and 
a 13 C nuclear spin (J = 1/2, g n = 1.4048) with an exter- 
nal magnetic field B, the second the zero- field splitting 
(ZFS) (D/h = 2.877 GHz, E/h = 7.7 MHz) of the NV 
ground state and the last the hyperfine interaction of the 
electron and nuclear spin (hyperfine coupling parameters 
A\\/h = 200 MHz, A±/h = 120 MHz). The Bohr magne- 
ton and the nuclear magneton are denoted as j3 e and /3 n . 
The eigenenergy spectrum was obtained by diagonalizing 
the Hamiltonian and confirmed with EasySpin [24]. 

Although the main axis of an NV center is along any of 
the four crystallographic (111) directions, for the (001) 
crystal plane there are only two magnetically inequiva- 
lent directions, labeled I and II in Fig. 1(b). Hence, ap- 
plying a magnetic field in an arbitrary direction within 
the (001) plane divides the whole ensemble of NV cen- 
ters into two magnetic subensembles due to the electron 
Zeeman effect. Consequently, each subensemble will be 
tuned into resonance with the CPW resonator at a dif- 
ferent magnetic field value (see Fig. 2(d)). For the CPW 
resonator, we chose a capacitively coupled A/2 niobium 
SC-resonator, fabricated on sapphire substrate. A me- 
ander geometry ensures a large overlap of the magnetic 



mode volume with the spin ensemble. At zero magnetic 
field and with the diamond on top, the resonator has a 
center frequency of w r /2i: — 2.701 GHz and a quality fac- 
tor Q sa 3200. Comparing these values with the empty 
and coupling limited resonator (uj r /2ir — 2.897 GHz and 
Q k, 4 x 10 4 ) shows that the diamond introduces consid- 
erable losses, possibly due to surface conductivity effects 
[25]. All subsequently shown measurements of the res- 
onator S21 -parameter are performed by using a vector 
network analyzer with a probe power corresponding to 
roughly 8 x 10 4 photons in the resonator. Lowering the 
power did not alter the signal, except for lowering the 
signal to noise ratio [26]. During all measurements the 
temperature was below 60 mK, ensuring a high degree of 
spin-polarization and a low thermal photon background 
in the frequency range of interest. 

A calculation similar to Verdu et al. [2] adapted for 
NV centers (see also Kubo et al. [14]) shows that the 
small mode volume of a CPW resonator leads to a cou- 
pling strength of go/2n « 12 Hz for a single NV center 
1 //m above the center conductor. This is not enough to 
exceed the losses k/27t sw 0.8 MHz of the resonator and 
the decoherence rate of the color center. However, by 
coupling a large number of spins N sw 10 12 collectively to 
the same electromagnetic mode of the resonator, an en- 
hancement of g co i /2tt s» yiVgo/27r ~ 10 MHz allows to 
enter the strong coupling regime. This can be modeled 
by a generalized Hamiltonian 

H = hw^a +^Y1 Wi i a ii + / 'X!- / " : ' T '/" + at<T iP' 
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where the index i runs over the different Zeeman 
subensembles I, II and j accounts for an inhomogeneous 
distribution of spin transition frequencies as well as 
coupling strengths gy. The operators er^ and cr^ de- 
scribe the spin raising and lowering operators, and a 
denote the creation and annihilation of a cavity photon. 
Although the presence of inhomogeneous broadening of 
the spin transitions has interesting consequences [27, 28] , 
the subsequent analysis assumes ujij = cjj. The effective 
coupling strength then reduces to g co i = \/J2i 9i ■ Fur- 
thermore, if the number of excitations meets n exc -C yN, 
the spin ensemble is well approximated by a system of 
coupled harmonic oscillators [29]. 

For strong coupling, an avoided crossing with a level- 
splitting of 2g co i is expected when the ensemble of spins 
is tuned into resonance with the cavity. Fig. 2(a) shows 
the transmission through the resonator as the magnetic 
field amplitude is varied for a field direction of <p = 45°. 
There, a large Rabi splitting of 2gij2-K — 18.5 MHz is 
observed, corresponding to a single Zeeman subensem- 
ble in resonance (uj r = ujL). The FWHM linewidths 
of the two peaks are found to be slightly different rj 
2ir x (4.4, 3.3) MHz at the avoided crossing. By applying 
a field in the [100] direction (ip = 0°), where the energy 
levels of the two subensembles I and II coincide, we can 
directly measure the -s/iV-enhancement. Since twice as 
many spins couple to the field mode, an enhancement 
compared to the case for ip — 45° of g co i ~ V%9i is ex- 
pected. Indeed we obtain an enhancement by a factor of 
1.42 ( 2g co [/2ir = 26.3 MHz). We also note an increase 
of the linewidth of the two peaks to s» 2ir x 7.3 MHz for 
ip = 0°. In case the magnetic field projections experi- 
enced by the two NV subensembles do not coincide (e.g. 
ip = +3°), two distinct avoided crossings are observed 
(see Fig. 2(b), (d)). 

Qualitatively one can understand the observed differ- 
ence in peak amplitude and peak width in Fig. 2(c) by the 
presence of higher lying NV states ( \+)i, |— )//). Includ- 
ing these states leads to a model of three (or more) cou- 
pled oscillators, where the peaks do not have equal height 
and width [30] . Furthermore inhomogeneous broadening, 
which is not contained in a model of coupled oscillators, 
influences the width and shape of the peaks [28]. 

Another interesting observation is the appearance of 
weak satellite lines around the main ESR line, visual- 
ized in Fig. 3(b) by the resonator transmission ampli- 
tude as a function of the magnetic field. As previously 
shown [23] , the hyperfine interaction of a nearest neigh- 
bor 13 C nuclear spin with the NV spin leads to an ESR 
doublet split by 130 MHz in zero field. However, for 
our geometry and a non-zero magnetic field, the pseudo- 
nuclear Zeeman effect [31] gives rise to two additional 
transitions, resulting in four allowed ESR lines (selection 
rules |Ams| = 1, |Amj| = 0). In Fig. 3(a) we show the 
position of the satellite lines as a function of the mag- 
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FIG. 3. (Color online) (a) Measured (dots) and simulated 
(solid lines) nearest neighbor 13 C satellite lines as a function of 
the magnetic field angle (colors denote the subensembles I and 
II). Also shown is the measured (center of the avoided cross- 
ing) and theoretical position of the main NV ESR (squares 
and dashed lines). The dotted line marks the angle at which 
the data in (b) is taken, (b) Amplitude of the lower (black) 
and higher (gray) frequency peak of the transmission signal 
obtained from Lorentzian fits as a function of the magnetic 
field. Satellite lines due to the hyperfine interaction of nearest 
neighbor 13 C nuclear spins with the electronic NV are visible 
and indicated by arrows. For p = 12.3° we see the four al- 
lowed ESR transitions of subensemble I (red arrow) and one 
of subensemble II (green arrow). 



netic field direction. These findings agree well with the 
simulation of the full spin Hamiltonian Eq. 1 when we 
assume hyperfine coupling parameters Au/h = 200 MHz 
and A±/h = 120 MHz [32]. Since the isotope 13 C only 
appears at its natural abundance of 1.1% in our samples, 
we estimate the coupling strength of these resonances to 
be weaker by a factor of y^3 x 1.1/100. For this reason 
we do not observe an avoided crossing for these lines but 
only a reduction of the transmission amplitude and an in- 
creased linewidth. Also note that some remote satellite 
lines observed are not predicted by the model. 

Additionally, we use a cavity readout scheme in the dis- 
persive regime of cavity QED, where A = \co r — uj±\ 3> 
\/~N g. Here the atomic ensemble can be seen as a re- 
fractive medium, lowering the resonator frequency by 
A c ~ Ng 2 /A when the ensemble is polarized in the 
ground state. If we add an additional MW tone u>d = oj± , 
a fraction of the NV ensemble is excited, leading to a 
positive shift which is plotted in figure 4. In Fig. 4(a) we 
scan the spectroscopy frequency 0Jd across the NV ESR 
lines and obtain the zero-field splitting parameters of the 
diamond sample by making use of the aforementioned 
dispersive interaction. We extract values for the parallel 
component of the zero- field splitting D/h — 2.877 GHz 
and the transversal component E/h = 7.7 MHz from 
the position of the minimum and the distance of the two 
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FIG. 4. (a) Dispersive resonator shift as a function of the 
spectroscopy frequency ujdf2^ (B = mT). This measure- 
ment allows to extract the zero-field splitting parameters 
D/h = 2.877 GHz and E/h = 7.7 MHz. (b) NV spin re- 
laxation measured via the resonator shift after the excita- 
tion frequency has been switched off (B = mT, ujd/2ir = 
2.885 GHz). Black dots refer to resonator frequency measure- 
ments after a given time after switch off (r = 44 ± 4 s), di- 
amonds mark measurements with a continuously probed res- 
onator (r = 45 ± 4 s) (see text for details). An exponential 
fit of the black data points is also displayed. 



peak maxima, respectively. As a fit function we used a 
simple asymmetric lineshape [33] . These values are in ac- 
cordance with ODMR measurements of the same sample 
and taking temperature effects into account [34]. 

We also present a relaxation measurement of the spin 
ensemble after excitation at w<j/27r = 2.885 GHz (see 
Fig. 4(b)). In this measurement we first excite the NV 
centers, then wait for a given time after the dispersive 
tone has been switched off and finally probe the cavity 
resonance. This sequence is repeated for different wait- 
ing times. The resonant cavity probe signal is switched 
on only in this last part of the measurement sequence. 
From an exponential fit we determine a relaxation time 
constant r = 44 ± 4 s. This value agrees well with previ- 
ous Tl measurements of NV centers in diamond at low 
temperature [35, 36]. In order to study the influence of 
the probe signal on the spin relaxation, we excite the 
NV centers as before and then probe the cavity contin- 
uously on resonance. For low probe powers ( ~ 5 x 10 4 
photons) the obtained time constant r = 45 ± 4 s agrees 
with the previously measured decay time. This shows 
that the measurement is non- destructive which is due to 
the fact that the dispersive measurement is a quantum 
non-demolition measurement at low probe power. How- 
ever, we note that for higher probe powers the observed 
time constant substantially decreases and the form of the 
decay deviates from a simple monoexponential decay. 

In conclusion we have shown the strong coupling of 
an NV spin ensemble to a CPW resonator and its vN- 



scaling with the number of emitters. We determined 
the spin-lattice relaxation time and ZFS parameters at 
60 mK using dispersive measurements, a powerful tech- 
nique in cavity QED. Furthermore the observed coupling 
to nuclear spin degrees of freedom in diamond may pro- 
vide an avenue towards extended storage times and quan- 
tum information processing with auxiliary nuclear spins. 
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